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Quantitative Pulmonary Imaging 



FIELD OF THE INVENTION 

The present invention provides a novel non-invasive method to quantitatively evaluate 
pulmonary ventilation and perfusion, combining the use of a hyperpolarized noble gas for 
5 ventilation MR images and arterial spin-tagging for perfusion MR imaging. Also provided are 
the algonthmic tools necessary to calculate the quantitative imaging of limg ventilation. 



REFEBENCE TO RELATED APPLICATIONS 

This application claims priority to 60/267,282 filed February 8, 2001, herein 
1 0 incorpo rated by reference in its entirety. 

GOVERNMENT INTERESTS 

This invention was supported in part by Grant Nos. K23 HL04486, RR02305, and ROl- 
HL-64741 from the U,S. National Institutes of Health. The Government may have certain 
1 5 rights r 1 this invention. 



BACKGROUND OF THE INVENTION 

The accurate measurement of pulmonary ventilation (V), perfusion (Q), and V/Q 
CO distribution is of great importance in medicine. V/Q relationships are used to diagnose and 

20 treat pulmonary disorders, including emphysema, asthma and pulmonary embolus. Established 
methods for anatomic study of gas exchange in the human lung include conventional planar 
ventilation and perfusion nuclear scanning (V/Q scans) and single photon emission-computed 
tomography (SPECT). However, existing techniques have well-documented, clinically 
signif cant hmitations. Moreover they are unfortunately limited in spatial resolution, and they 
25 involve the mandatory exposure of both the patient and the practitioner to radioactive agents. 
They also require prolonged breathing from an increased resistance closed system through a 
tight-fitting facemask, an effort that some patients with advanced emphysema are unable to 
tolerate. 

Normally, pulmonary artery blood flows through capillaries that are immediately 
30 proximate to alveoli. For gas exchange to occur efficiently, alveolar ventilation and pulmonary 
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blood flow are proportionally matched. When V/Q mismatch occurs, the terms "shunting" or 
"dead space" are used to describe the abnormalities in V or Q, respectively. Shxmting occurs 
when unventilated or collapsed alveoU cannot participate in gas exchange (/.e., unloading 
carbon dioxide and absorbing oxygen), so that there is diversion in the blood flow. Likewise, 
5 without blood flow to an alveolus, wasted ventilation or dead space occurs. While small 

amounts of shunting (e.g., in the bronchial vessels) and dead space (e.g., in the nasopharyx) are 
normal, a significant number of pathological pulmonary conditions are characterized by V/Q 
mismatches. 

Pulmonary functional imaging methods include qualitative tests that provide valuable 
10 clinical information, but each currently available test has a limitation. The radionucUde V/Q 
scan has inherently poor spatial relations, lacks the abihty to image cross-sectional views, and 
deposits radiolabeled aerosols in the central airways of the patient, causuig errors. SPECT has 
the highest resolution of all of the known methods, but accurate detection of pulmonary 
13 disorders still demands even higher resolution than is currently provided. Global pulmonary 

i2 15 function may be determined by the multiple inert gas elimination technique (MIGET), which 
4^ provides a quantitative measurement of pulmonary V/Q ratios based upon a mathematical 

subdivision of the lung into fifty hypothetical segments. It is the gold-standard physiological 
f test of lung V/Q relationships, but even MIGET has characteristic and problematic limitations, 

f J which have prevented its widespread clinical appUcation. For instance, large blood sample 

H 20 volumes are necessary for the analysis. 

Ventilation scanning using non-radioactive xenon gas and ultrafast computed 
tomography has been described (Murphy et al. Chest 96:799-804 (1989), but requires exposure 
to radiation in the form of x-rays. Low signal intensity and the effects of respiratory motion 
have likewise limited conventional proton MR imaging of human airways and the pulmonary 
25 parenchyma (Bergin et al, Radiology \19\111-1%\ (1991); Bergin et al, Radiology 180: 845- 
848 (1991); Bergin et al, J, Thorac. Imag. 8:12-17 (1993); Alsop et al, Magn. Reson. Med. 33: 
678-682(1995)). 

The use of hyperpolarized noble gases, such as ^He, has been demonstrated to be useful 
in the imaging of gas distribution (ventilation) m the human lung (Middleton et aL, Mag. Res, 
30 Med. 33:271-275 (1995); MacFall et al, Rad 200:553-55% (1996); Kauczor et al, Rad 

201:564-568 (1996);. Kauczor et at, J. Mag Res. Imag. 7:538-543 (1997); Roberts et a/., Mag 
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Res. Med. 44(3):379-382 (2000); Black et ah, Rad 199(3):867-870 (1996); de Lange et al. Rod 
2I0(3):851-857 (1999); Altes a/., J. Mag. Res. Imag. 13(3):378-384 (2001); Salerno etal, 
Mag. Res. Med. 46:667-677 (2001); Gierada et al, NMR Biomed. 13(4):176-181 (2000)). 
When hyperpolarized heliiim-3 (HP-^He) has been used in the lung, regional oxygen partial 
5 pressure (O2) has been calculated from the intensity of the HP- He magnetic resonance (MR) 
images, but it has not been previously applied to V/Q measurements. 

Thus, the ideal radiologic method for assessing gas exchange would quantify the 
comparison of V/Q in small volumetric zones of the lung. It would be safe, non-invasive, 
rapid, easily applied, reproducible and tolerable to patients with advanced lung disease. To 

10 date, no available technique met all of these requirements. 

Arterial blood perfusion imaging can be obtained by noninvasive magnetic labeling of 
blood using the arterial spin-tagging (AST) technique (see, e.g., U.S. Patent No. 5,402,785, 
Leigh and Alsop). This method is based upon the steady-state magnetic labeling, or "spin- 
tagging" of arterial water protons (Detre et al,, NMR Biomed, 7:75-82 (1994); Roberts et al, 

15 Proa Natl Acad. ScL USA 91: 33-37 (1994); Roberts et al, Proc. Fifth Scientific Meeting and 
Exhibition of ISMRM 3:1764 (1997); Roberts etal, i?aJ212(3):890-895 (1999); Detre etaL, 
Magn, Res. Med. 22:37-45 (1992); Roberts et al, 1 Magn. Reson. Imag. 14(2):175"180 
(2001)). 

The combined application of these two non-invasive functional MR imaging techniques 
20 would provide images of the patterns of air flow in the lungs of human or animal subjects, and 
aid in V/Q mapping of the lung. Prior to the present invention, however, although qualitative 
methods are known, no noninvasive quantitative method for imaging lung ventilation has been 
used or developed. Moreover V/Q by magnetic resonance imaging (MRI) would offer a new 
ability to uniquely describe normal lungs, to diagnose pulmonary disease, to follow the course 
25 of pulmonary diseases, and to assess therapeutic response in terms of localized V/Q (regional 
V/Q). Accordingly, a need has remained in the art until the present invention for a quantitative 
method of rapidly providing accurate V/Q images at very high resolution, and an invention 
meeting these needs would provide valuable insight into the pathogenesis and pathophysiology 
of pulmonary V/Q abnormalities and of the normal lung. 
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SUMMARY OF THE INVENTION 

The present invention offers novel quantitative methods and tools for computing three- 
dimensional quantitative images of human lung ventilation. Also provided are quantitative 
methods for very high resolution, three-dimensional combined imaging of pubnonary 
5 ventilation (V) and perfusion (Q), synchronously (V/Q), wherein the methods comprise 

delivering a predetermined volume of hyperpolarized noble gas into the conducting airways in 
each ventilated region of the pulmonary system and collecting local magnetic resonance image 
data therefrom to produce a V/Q image data set. 

Using an inhaled noble gas, preferably laser-hyperpolarized helium-3 (^He) gas (HP- 
10 ^He), three-dimensional MRI images are rapidly obtained. The HP-^He acts as a contrast agent, 
thereby providing unparalleled, high resolution images defining the air spaces. 

The methods comprise three primary components, including the algorithmic tools 
^ therefor, each of which is provided in the present invention. First, the methods rely upon 

O combined MRI ventilation/perfusion (V/Q) imaging techniques, which are provided and 

|i 15 applied for the first time synchronously in humans in the present invention. The quantification 
■^'r- of regional V/Q is extremely important for diagnosis and evaluation of pulmonary diseases, 

:=ip A second provided component is the three-dimensional quantitative imaging of absolute 

f 5 lung perfusion (Q) and collecting local magnetic resonance image data therefrom to produce an 

2 absolute lung perfusion image data. In the present invention quantitative imaging of absolute 

tg 20 lung perfusion is accomplished by any known method, including gadolinium-based methods. 

However, in a preferred embodiment of the invention, Q is computed using the arterial spin- 
tagging (AST) method of perfusion MRI, 

A third provided component for the methods of the present invention is the co- 
registration algorithm to spatially co-register the two image data sets, (HP-^He MRI image of 
25 V/Q and MR imaging of quantitative perfusion (Q) in the lung). The images of V/Q and 
perfusion (Q) are generally obtained at separate times and during different phases of the 
respiratory cycle. For this reason, a spatial co-registration algorithm of the present invention 
must be used to align the image data prior to calculation of ventilation. 

From the data acquired in the combined data sets and their spatial co-registration 
30 resulting from the measurements and calculations of the present invention, absolute ventilation 
(V) is computed by multiplying the ventilation/perfusion ratio (V/Q) by perfusion (Q) for each 
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point in the lung image, thereby producing three-dimensional quantitative images of Ixmg 
ventilation. 

The methods of the present invention accomplish the quantitative calculation of 
absolute lung ventilation by dividing the lung image into as many distinct voxels as the 
5 imaging resolution permits. Regional PO2 (partial pressure of oxygen gas) and PCO2 (partial 
pressure of carbon dioxide gas) measurements resulting from the HP-^He MR images enable 
the calculation of regional V/Q distribution, and in particular alveolar V/Q distribution, by 
MRI. Then, once the ventilation/perfusion ratio (V/Q) is determined, it is multipUed by the 
perfusion value (Q) for each point in the lung image. 
10 Further provided in the present invention are the very high resolution, three-dimensional 

images of pulmonary ventilation (V) and perfusion (Q) in the lung, produced by the foregoing 
methods using hyperpolarized HP-^He gas. 
. In addition, systems are provided for producing the very high resolution, three- 

dimensional images of pulmonary ventilation (V) in the lung in accordance with claim 8, 
15 comprising: means for collecting and processing magnetic resonance imaging data from the 
hyperpolarized H^He gas-infused lung, conducting airways and ventilated regions of the lung; 
means for dividing the lung images into as many distinct voxels as imaging resolution permits; 
means for synchronously imaging puhnonary ventilation (V) and perfusion (Q), and calculating 
V/Q therefrom; means for imaging of absolute lung perfusion and calculating Q therefrom; 
: 3 20 means for spatially co-registering the HP-^He MR image V/Q data set with data set Q for 
[ t quantitative lung perfusion; and means for computing absolute ventilation (V) from the 
= "nri ventilation/perfusion ratio (V/Q) and perftision (Q). 

A device is also provided for post-acquisition processing of MRI-acquired high 
resolution, three-dimensional images of pulmonary ventilation (V) in the lung, comprising: 
25 computer-readable signal-bearing medium; means for collecting and processing magnetic 

resonance imaging data from a hyperpolarized H^He gas-infused lung, conducting airways and 
ventilated regions of the lung; means for dividing the lung images into as many distinct voxels 
as imaging resolution permits; means for synchronously imaging pulmonary ventilation (V) 
and perfusion (Q), and calculating V/Q therefrom; means for imaging of absolute lung 
30 perfusion and calculating Q therefrom; means for spatially co-registering the HP-^He MR 
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image V/Q data set with data set Q for quantitative lung perfusion; and means for computing 
absolute ventilation (V) from the ventilation/perfusion ratio (V/Q) and perfusion (Q). 

Additional objects, advantages and novel features of the invention will be set forth in 
part in the description, examples and figures which follow, and in part will become apparent to 
5 those skilled in the art on examination of the following, or may be learned by practice of the 
invention. 



DESCRIPTION OF THE DRAWINGS 

The foregoing summary, as well as the following detailed description of the invention, 
10 will be better understood when read in conjunction with the appended drawings. For the 

purpose of illustrating the invention, there are shown in the drawings, certain embodiment(s), 
which are presently preferred. It should be understood, however, that the invention is not 
1^ limited to the precise arrangements and instrumentalities shown. 

:S FIG. 1 depicts the three-dimensional compact, gradient-echo pulse sequence used to 

15 acquire perfusion image data showing the radiofrequency (RF) and gradient (Gx, Gy, and Gz) 
waveforms. An arterial labeling RF pulse was applied in the presence of a gradient to provide 
continuous magnetic labeling of pulmonary arterial blood flow. This pulse had an amplitude 
(Bi) of 24 mgauss and a duration (Ti) of 40 mseconds corresponding to a duty cycle of 80%. It 
was appUed during the time between the crusher/refocusing gradient pulses (C) and the 
20 beginning of next excitation. The flip angle of the excitation pulse, a, was 30°. The time axis 
is not to scale. 

FIGs. 2 A-2C depict the lung of a healtiiy male nonsmoker. FIG. 2A depicts a sagittal 
proton image of the lung. The corresponding images were depicted showing ^He distribution 
(ventilation) in FIG. 2B, and perfusion via the AST technique in FIG. 2C. 
25 FIGs. 3 A-3D depict corresponding proton (FIG, 3A) and ^He (FIG. 3B) coronal images 

of a patient with a history of emphysema. The patient had received a left lung transplant, and 
the ^He distribution in the transplanted left lung allograft was uniform. The *^He ventilation 
images show large defects in ^He distribution in the right (native) lung. Corresponding AST 
sagittal images perfiision images from the right native lung (FIG. 3C) show large defects in 
30 signal, while images of the left lung allograft (FIG. 3D) show bright signals from normal blood 
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perfusion. Multiple slices may be obtained. Three images of the chest are shown for each view 
to cover a broad spatial range. 

FIG. 4 depicts sagittal images of pulmonary perfusion obtained using the AST method 
in a patient with a large puhnonary embolus, documented by computer tomography (CT) in the 
5 right lower lobe puhnonary artery. The marked decreased perfusion signal in the lower lobe is 
shown with arrows. As in FIG. 3, multiple slices may be obtained of the lung. In this case four 
images are shown, but the selected number of images is not relevant. 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
10 The present invention comprises imaging and quantitative measurement of lung 

ventilation, particularly in a human lung. Methods for quantitative imaging of lung ventilation, 
and algorithmic tools therefor, comprise three primary components, each of which comprises 
an element of the present invention, and which will be described in greater detail below. 
First, the methods rely upon combined MRI ventilation/perfixsion (V/Q) imaging 
15 techniques, which are provided and applied in combined synchronous form in a patient for the 
first time in the present invention. Using an inhaled noble gas, preferably laser-hyperpolarized 
helium-3 (^He) gas (HP-^He), three-dimensional MRI images are rapidly obtained of V/Q in the 
lung. 

A second component comprises the three-dimensional quantitative imaging of absolute 
20 lung perfusion (Q), such as by the gadolinium based methods, or any known method, so long as 
V and Q imaging are combined. In a preferred embodiment, perfusion is measured using the 
arterial spin-tagging (AST) method of perfusion MRI, as has been previously described in, e.g., 
U.S. Patent No. 5,402,785 (herein incorporated by reference in its entirety). 

A third component comprises the co-registration algorithm to spatially co-register the 
25 two image data sets (the HP-^He MRI images of V/Q and the quantitative perfusion (Q) in the 
lung, preferably by AST MRI). Since the images of V/Q and perfusion (Q) are generally 
obtained at separate times and during different phases of the respiratory cycle, a spatial co- 
registration algorithm is needed to align the image data prior to calculation of ventilation. 
Finally, absolute ventilation (V) is computed by multiplying the ventilation/perfusion ratio 
30 (V/Q) by perfusion (Q) for each point in the lung image, thereby producing three-dimensional 
quantitative images of limg ventilation (in units of ml*g"^*min"^). 
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Preferred embodiments of the present invention produce high-resolution, three- 
dimensional images of puhnonary ventilation and perfusion (V/Q) in human and animal 
subjects using a noninvasive and nonradioactive MRI technique. The novel tools and methods 
for its use to evaluate pulmonary V/Q combine the use of a hyperpolarized noble gas for 
5 ventilation images and AST for perfusion imaging. 

The present invention offers several advantages over existing methods for studying V/Q 
relationships in the human lung. The technique does not involve exposure to ionizing radiation. 
Thus, it avoids the risks associated with radioactivity, for both the patient and the clinician. It 
is completely noninvasive and, therefore, avoids the costs and risks associated with prior art 
10 methods involving gadolinium administration. However, use of the AST methods will not 
degrade subsequent gadolinium-based studies, such as magnetic resonance angiography. 
Additionally, multiple repeated images may be acquired by the present method throughout the 
day with no degradation of image quality. In contrast to radionuclide methods, the ability to 

Q 

C3 repeat the imaging process within a short time period in the same subject permits the 

^ f 15 measurement of dynamic responses to various stimuh. 

:? Other potential advantages include the fact that the ventilation portion of the study 

,C requires only a short breath-hold, while the AST perfusion technique does not require any 

patient breath-holding. This is a major advantage because many patients with pulmonary 
f U disease are simply unable to cooperate in breafh-holding maneuvers, and it also makes the 

20 process useful for animal subjects. 

This is the first time that these ventilation and perfusion techniques have been used 
simultaneously in patients to obtain functional magnetic resonance V/Q images. As such, the 
present invention offers great promise for evaluating apical perfusion in patients before and 
after lung volume reduction surgery Kotloff aL, Chest 120(5):1609-1615 (2001), or 
25 potentially for detecting puhnonary emboli. These techniques also are useful in the evaluation 
of patients before and after the use of bronchodilators or vasoactive substances. Moreover, this 
form of functional pulmonary imaging paves the way for noninvasive in vivo V/Q mapping of 
the human lung. 

Imaging of lung ventilation/perfusion (V/Q) ratio using laser-polarized Helium-3 gas. 
30 In the preferred embodiment of the method of the present invention, V and Q imaging is 

performed sequentially. The helium imaging is performed during a breath-hold to obtain gas 
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distribution (ventilation) images, then the helium is exhaled. Perfusion imaging is then 
performed during normal breathing, and the AST method allows for the perfusion information 

to be obtained while the patient is breathing normally. Perfusion data can also be obtained 
using gadolinium based methods, or any known method, so long as V and Q imaging are 
combined. 

In practice, the patient's lungs are ventilated, and imaging is performed using 
hyperpolarized ^He (HP-^He). The polarized gas was inhaled by the patient while lying in a 
supine position in the MRI magnet device. The subjects inhaled approximately 750 cc of the 
HP- He gas to near total lung capacity, and afler adequate mixing, the lung was imaged with a 
homogeneous coil for both signal transmission and reception, as described in the examples that 
follow. 

Because time of repetition (TR) is irrelevant to polarized gas imaging, a real-time series 
of images is obtainable. This means that V/Q by MRI will be able to uniquely describe any 
focal or widespread pulmonary disease, and to define the response to therapy in terms of 
regional V/Q. For example, the MRI techniques associated with the present invention are 
particularly suited to describe alleviation of V/Q mismatch after prone positioning in adult 
respiratory distress syndrome (ARDS). To the medical community, functional lung studies and 
quantification of regional V/Q relationships have far-reaching clinical significance, offering 
unprecedented means to diagnose and follow the progression of puhnonary disease, as well as 
to regionally track the patient's response to therapy. 

In a preferred method, the collected lung images are divided into as many distinct 
voxels as the imaging resolution permits, and the starting point for calculating local lung 
volume is the average signal intensity in each voxel divided by the signal intensity in the 
trachea. Neglected from the calculation are: (1) any signal decay due to oxygen or wall 
interactions; (2) any loss of heUum from the gas spaces due to absorption by tissue or blood; 
and (3) any contribution to the signal from the tissue space within the voxel. The average 
concentration of HP-^He in each voxel is calculated in terms of the amount of HP-^He in each 
voxel divided by the volume of the voxel. The amount of HP-^He in the voxel is simply 
calculated in terms of the concentration of HP-^He in the gas space of the voxel multiplied by 
the volume of the gas space in the voxel. Assuming complete equilibration of the HP-^He 
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throughout the lung, the concentration of HP-^He in the gas space of the voxel is exactly equal 

to the concentration of HP-^He in the trachea. 

MRI imaging with HP-^He is ideal for the three requirements of the preferred 

embodiment of the method of the present invention - high resolution, high signal intensity, and 
5 rapid time resolution. Commercially-available MRI systems, such as those manufactured by 

GE and Siemens, and improvements thereof are well known in the art, and may be used and 

adapted for these methods. 

The subject of the present method is a mammal. Preferably, the manmialian subject is 

human. The lung(s) of such subject may be either normal for mapping purposes, or injured or 
10 diseased for therapeutic or diagnostic purposes. Both halves of the subject's lung need not be 

in the same condition, and may be compared against each other. Alternatively, a normal lung 

may also be used for comparison purposes. 
H Regional PO2 (partial pressure of oxygen gas measurements) computed fh>m HP-^He 

MR rniages enable the calculation of regional V/Q distribution by MRI. Magnetic image data 
■'^ 15 are transferred by a high-speed network to a workstation computer (Silicon Graphics, Inc.), and 
£ image analysis is performed using IDL (Research Systems, Inc. Boulder, CO). In the ^He 

■ S images, the normaUzed first breath ^He untake, U\e, was computed by placing a ROI (Region 

of Interest) interactively around the trachea. The mean signal intensity in this ROI is measured 

and the magnitude images are scaled by &e resulting value, 
y 20 The decay of the HP-^He polarized signal in a series of lung images can be modeled by 

C3 taking into account three combined effects. Specifically, the loss of the Helium-3 

magnetization ("signal") results from three effects: (i) radiofrequency excitation pulses ('B 1', 

coil-dependent), (ii) the presence of molecular oxygen (O2), and (iii) interactions with the 

alveolar wall. 

25 Previous work by the Deninger group in Mainz, Germany (Deninger et aL, J. Magn. 

Reson. 141(2):207-216 (1999); Deninger et al, NMR Biomed. 13(4): 194-201 (2000)) has 
shown that the wall interactions (effect (iii), above) are negligible. Hence, the reduction in 
signal intensity is mostly due to: (i) the RF excitation on the nuclear spin polarization of HP- 
^He in the acquisition volume, and (ii) decay in signal intensity due to relaxation between 

30 subsequent images, which is a function of the local oxygen pressure. The ratio of signal 

intensity in two consecutive images in a given region of interest depends not only on factors (i) 
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and (ii), as outlined above, but also on the time interval (At) between subsequent acquisitions. 
By incorporating the combined effects of these factors on spin relaxation mechanisms, it is 
possible to model the dynamics of the signal in a series of lung images (Deninger et aL, 1999), 
and thereby obtain the following equation: 



in which a is the flip angle of the RF pulse, q is the ratio of signal intensity of two consecutive 
images of S^^^ and 5„ , p^^ , and r is the number of phase encoding steps in each image. Since 
the p^^ is constant within an imaging experiment, obtaining different measurements at constant 
Pq^ and different flip angle a leads to the following equation: 



Thus, with known a, it is possible to calculate the regional [O2] of oxygen in the lung. 
This method, to our knowledge, was first described and published by the Deninger group 
(Deninger et al, 1999; Deninger et al, 2000). 

Next, the values for the V/Q ratio from the PO2 measurement are computed. The 
average Va/Q ratio (where Va is alveolar ventilation) is calculated in each voxel from the 
measurements of average PO2 per voxel, using well-described equations of steady state 
exchange of O2, CO2, and N2 to the local Va/Q ratio. This is a variation of established 
techniques that calculate Pa02 (alveolar PO2) from the Va/Q ratio. Calculation of alveolar PO2 
(PAO2) from an assumed Va/Q ratio is known as the "forward problem," Detailed techniques 
for performing this calculation have successftilly predicted overall pulmonary gas exchange 
from measured Va/Q distributions. Identical techniques are used to calculate Va/Q ratios from 
measured PO2. 

First the numerical methods will be described for solving the forward problem, then the 
necessary modifications will be described for calculating Va/Q from the local PO2 
measurements. Methods for calculating PO2 and PCO2 (partial pressure of carbon dioxide gas) 
in a homogeneous lung unit as a fimction of the Va/Q ratio are well estabhshed, and are 
commonly used to predict overall physiologic gas exchange for a given Va/Q ratio distribution. 
In the forward problem, it is assumed that Va/Q is known, as a result, lung inputs, such as 
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inspired oxygen fraction (F1O2), mixed venous PO2 (Pvo2), and mixed venous PCO2 (Pvco2). 
From this assumed YJQ ratio, the object of the forward problem is to calculate the unknown 
alveolar gas tensions Pa02 and PaC02. These calculations start with steady state mass balances 
for gas exchange. The result is four equations in four unknowns: (1) Pa02, (2) PaC02, (3) Vi/Q 
5 (ratio of inspired ventilation (Vi may not equal Va) / local perfusion (Q)), and (4) PaN2 

(alveolar nitrogen partial pressure). By elimination these are reduced to two equations in two 
unknowns, Pa02 and PaC02. In particular, the equations relating to blood partial pressure to 
blood content for O2 and CO2 are highly nonlinear. Therefore, the two equations are expressed 
as two imphcit functions of the two unknowns, and the system is solved iteratively using the 
10 secant method (see, e.g., Numerical Recipes in C , WH Press (Teukolsky, Vetterling & 
Flannery, ed.), Cambridge University Press, 1992). 

Va/Q was obtained from measured PO2. The starting point is again four equations in 
four unknowns that result from the mass balances for steady state gas exchange. These four 
H equations are reduced to two equations in two unknowns Va/Q and PaC02. In the inverse 
H 15 problem, the system is further reduced to one equation and one unknown. This equation again 
V is expressed as an implicit function of the unknown PaC02, and the numerical solution is 
^ obtained by the secant method. Steady state mass balance on oxygen exchange results in 

Equation 3. 



wherein Pi02 is inspired oxygen partial pressure, k is a constant to convert gas volumes, Cc'02 is 
end-capillary O2 content, and Cvo2 is mixed venous O2 content. Cvo2 is calculated from 
numerical O2 and CO2 dissociation curves (see e.g., Olszowka et al. In: Puhnonary Gas 
25 Exchange (1st ed.), JB West (ed.), Academic Press, New York, New York, 1980, p. 263-306. 
The measured mixed venous PO2, PCO2, and pH. Cc'02 is calculated from these curves, but the 
unknown PaC02 is an imphcit variable. Steady state mass balance on CO2 exchange results in 
Equation 4. 



20 



((Vi/Q) X (P1O2)) - ((Va/Q) X (Pa02)) = k(Cc.02) - (Cvo2) 



(Equation 3) 



30 



((Va/Q) X (PACO2)) = k(Cvco2) - (Cc'co2) 



(Equation 4) 
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wherein Cc co2 is end-capillary CO2 content, and Cvco2 is calculated from numerical O2 and 
CO2 dissociation curves (see, e.g., Olszowka et al, 1980, supra). For N2, a mass balance leads 
to Equation 5. 

5 ((Vi/Q) X (P1N2)) - ((Va/Q) X (PaN2)) = h^i (PaN2) - (PvN2) (Equation 5) 

wherein P1N2 and PaN2 are inspired and alveolar nitrogen partial pressures, respectively, A,n2 is 
the blood gas partition coefficient for N2, and Pvn2 is the mixed venous nitrogen partial 
pressure (calculated from a balance of mixed venous blood gas data). The requirement for all 
10 gas partial pressures to sum to barometric pressure leads to Equation 6. 

(PAO2) + (PACO2) + (PaN2) = Pb - Ph20 (Equation 6) 

wherein Pb is barometric pressure and Pmo is water vapor pressure. Equations 5 and 6 are 
15 linear, and PaN2 and Vi/Q are eliminated by substitution, resulting in 2 equations, and in the 
two unknowns Va/Q and PaC02. One equation is Equation 4. The second equation, after 
substituting appropriate expressions for F1O2 in place of P1O2 is Equation 7. 

20 (Va/Q)x[PA-Fx0.xP.CO.-PaO.)] ^kCCcoa- Cvco^) (EquationT) 

(1 - F1O2) + A,N2 (Pb - Ph20 - Pa02 - PaC02 - Pvn2) 

In the forward problem, these two equations include implicit functions of both 
unknowns, Pa02 - PaC02. In the inverse problem, the additional simpUfication is that in both 
25 equations, fimctional dependence on Va/Q is explicit. Therefore, the system is simplified by 
eliminating Va/Q to give a single equation in one unknown, which is rearranged in the form 
F(PaCO2) = 0: 

[k(Cvco2 - Cc.co2) / PACO2] [P1O2 - F1O2 X PaC02 - Pa02)] / (1 - F1O2) + Xn2 (Pb - Ph20 - 
30 Pa02 - PaC02 - Pvn2) - k(Cc'02 - Cvo2) = 0 (Equation 8) 

This equation contains two implicit and nonlinear fimctions of PaC02 (Cc'co2 and Co'02)- 
It cannot be solved directly, but is solved iteratively by the secant method. The resulting 
PaC02 is uiserted into either Equation (2) or Equation (5) for direct calculation of Va/Q. 
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Imaging of absolute lung perfusion (Q) using arterial spin-tagging method of perfusion 

MRI. 

The "arterial spin-tagging" (AST) methods for quantitative perfusion MRI are known 
(see, e.g,, U.S. Patent No. 5,402,785, Leigh and Alsop), and are exemplified in the preferred 
5 embodiments of the invention. However, perfusion data is also provided by gadolinium based 
methods, or by any recognized method presently known or yet to be developed, so long as the 
V and Q imaging methods are combined. A brief summary of the AST approach, as modified 
to perform lung perfusion imaging follows. 

In the AST perfusion method, continuous magnetic labeling of the arterial water results 
10 in a small signal change that is a function of perfusion, arterial-capillary transit time, and 

regional blood volume. As shown by Zhang et ah {Magn. Reson. Med. 25:362-371 (1992)), a 
simplified model can be used to analyze the AST method. In this model, the Bloch equation 
(Detre et al., Magn. Reson, Med, 22:\-9 (1991)) is modified to include terms for perfusion, /, 
and exchange with a macromolecular pool and is given by: 



i = Vl ' -^forH -^Key^m + /^a - /^^ (Equation9) 



dM 

Tr=— TT^ 



in which: 



: :| / = tissue-specific perfusion in ml*g"l*sec""l, 

C3 T 1 J = spin-lattice relaxation time of tissue water, 

W 20 = longitudinal magnetization of water per gram of tissue, 

in M? = equilibrium value of Mt, 

s Li 

= longitudinal magnetization of water per ml of arterial blood. 
My = longitudinal magnetization of water per ml of venous blood, 

= longitudinal magnetization of macromolecules per gram (g) of tissue, 
25 kfoj., k^gy = magnetization transfer rates between free tissue water and 

macromolecular water. 

The tissue magnetization is measured under steady-state conditions of arterial labeling 
(Mj " ) and under control conditions in which the blood flow is not labeled (M^^°"), Solving for 
30 the observed alteration in magnetization in terms of perfusion, /, leads to Equation 10: 



dM 
M 



= 7'lap,exp 



la 



(Equation 10) 
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where 5M/M is the observed fractional change in magnetization induced by arterial labeling, a 
is the degree of magnetic labeling of the inflowing blood, xi is the arterial-capillary transit time, 
rib is the longitudinal relaxation time of the arterial blood, X is the tissue:blood partition 
coefficient, and Tl^^^^ the "apparent 71", is defined by: 



1 



1 / 



Tl 



(Equation 1 1) 



in which Tit is the longitudinal relaxation time of the tissue. This model is then modified to 
include a term for the blood volume (Roberts, "Magnetic Resonance Imaging of Perfusion in 
Humans using Spin-Tagging of Arterial Water," Thesis, 1994, University of Pennsylvania, pp 
58-60), yielding: 



dM 
M 



TL 



(Equation 12) 



in which v is the fi-actional blood volume and T2 is the transit time required for blood to reach 
the blood volume compartment. Equation 12 shows that the signal changes observed upon 
arterial labeling are a function of the tissue-specific perfusion, the fi^actional blood volume, the 
relaxation characteristics, and the transit times. 

A variant of the AST method has recently been described which is relatively insensitive 
to variations in the transit time and regional blood volume. In this approach, continuous 
labeling is followed by a delay time, co, during which unlabeled blood is allowed to "fill" the 
intravascular compartment at which point imaging occurs. Alsop et al {Radiology 208(2):410- 
416 (1998) and/. Cerebral Blood Flow & Metabol 16(6): 1236-1249 (1996)) have shown that 
in this approach perfusion may be calculated as: 



dM 
M 



exp 



la\ 



f 



(Equation 13) 
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in which w is the acquisition delay. This modification, termed *transit-time insensitive AST' 
(TTI-AST) avoids the problems associated with measurement of the transit time and regional 
blood volume. 

Spatial co-registration of the image data. 
5 The images of V/Q and perfusion (Q) are typically obtained at separate times and 

during different phases of the respiratory cycle. For this reason, a "spatial co-registration" 
algorithm must be used to align the image data prior to calculation of ventilation. 

Methods to solve the co-registration problem are known in the art. For example, one 
such method is based on the use of recognizable features within each data set, or "fiducials," to 
10 create a lattice of points. These points are then spatially-ahgned using a "mesh-warping 
algorithm." 

Specifically, the approach involves first globally aligning both MR images (in the 
following description, referred to as A and B). This comprises the following steps: (Step (a)) 
interpolate A and/or B so that they have the same voxel/pixel size. Let Al and Bl be the 
15 resulting images. (Step (b)) Roughly ahgn Al with Bl so that the subsequent elastic matching 
step will not have to worry about global aUgnment. Several methods are known in the art that 
satisfy this step and would be recognized by a practitioner in the field of image computing and 
calculation. 

The method used is based on the global mutual information of the joint histogram of Al 
20 and Bl . Such methods have been shown to yield very accurate and robust rigid registration, 
even when the relationship between the intensity patterns in the two images is very complex. 
Let A2 be the resulting image that is aligned in this global fashion with Bl, 

(Step (c)) This step performs elastic matching by local deformation. A2 is first divided 
into several sub-regions. The size of the sub-regions is and should be related to the degree of 
25 deformation that is expected between A2 and BL Among these sub-regions, those containing 
the recognizable "fiducials" (structural information) were selected, e.g., lung boundaries or 
pulmonary vessels (which are identified either interactively or by segmentation). 

Elastic registration is performed by optimizing the local image intensity correspondence 
of each selected sub-region as a function of its translation over B 1 . The required probabilities 
30 are obtained firom the joint histogram of A2 and B2. 
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(Step (d)) Finally the previous step (Step (c)) yields an elastic transformation field. 
This field comprises a vector associated with every voxel v in A2 that indicates how v must be 
moved to obtain the corresponding voxel in Bl. This information is then utilized to 
'\mdeform" A2 and create an image, A2 = B2, that matches (voxel-by-voxel) image Bl. 

Finally, after performing the measurements and calculations comprising the steps 
above, the "three-dimensional quantitative images of lung ventilation" (in units ml*g'^*min"^) 
are computed. 

The present invention is further described in the following examples. These examples 
are not to be construed as limiting the scope of the appended claims. 



EXAMPLES 

To confirm the effectiveness of the disclosed method, experiments were performed after 
informed consent was obtained from the patients using a medically approved protocol. Six 
patients were imaged. Three normal subjects, two patients with advanced emphysema, and one 

15 woman with a pulmonary embolus (PE), documented by history and V/Q scan, were selected to 
undergo magnetic resonance imaging. The two emphysema patients had each previously 
undergone single orthotopic lung transplantation (6 months and 5 years prior to imaging). 
Neither transplant patient had a baseline oxygen requirement, and both were firee of respiratory 
symptoms at the time of imaging. Single lung transplant patients were chosen so that a direct 

20 comparison could be made between a severely obstructed and a relatively normal lung in the 
same patient. In all of the patients, chest radiographs showed no infiltrates or other active 
disease. The patient with PE was receiving intravenous heparin at the time of imaging. 
Ventilation imaging was not performed in this patient as the diagnosis of PE was confirmed 
with the AST perfusion scan alone. 

25 In all subjects there were no apparent adverse effects firom inhalation of the 

hyperpolarized helium gas. Blood pressure and oxygen satiu'ation was continuously measured 
in the patients throughout the experiments. Oxygen saturations remained at all times greater 
than 96%. 

Magnetic resonance data were acquired at the MR imaging (MRI) center of the Hospital 
30 of the University of Pennsylvania using a superconducting 1.5 Tesla imaging system equipped 
with high-powered gradients (Echospeed Signa, 5,6 revision software, GE Medical Systems, 
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Milwaukee, WI). Proton imaging was performed at 63.8 MHz using a standard birdcage body 
coil for both signal transmission and reception. Helium imaging was performed at 48.65 MHz 
using a custom-built, broadband 26 cm diameter octagonal transmit/receive surface coil placed 
on the anterior thorax. 

5 Pulmonary ventilation imaging was perfonned using hyperpolarized ^He (HP-^He). 

The hyperpolarized gas was prepared through spin-exchange collisions with optically-pxmiped 
rubidium atoms. Helium was prepared by the metastability-exchange technique for one subject 
(Colegrove et al, Phys. Rev, 132:2561-2572 (1963); Gentile et al, Phys. Rev, A. 47:456-467 
(1993); Happer et al, Phys. Rev, A, 29: 3092-31 10 (1984); Ebert et aL, Lancet 347:1297-1299 
10 (1996); Saam, Nat Med. 2:358-359 (1996)). Polarization levels of 5%-15% were achieved. 
After polarization, the gas was transferred to a plastic balloon that was then given to the subject 
lying in a supine position in the MRI magnet device. The subjects inhaled approximately 750 
cc of the HP-^He gas to near total lung capacity and breathing was held for approximately 15- 
g 30 seconds to permit imaging of the ventilation while the HP-^He gas was held in the lungs and 
\! 15 pulmonary system» Breathing was resumed for the perfusion analysis. Once the breath-hold 
" Z has been exhaled, no further ventilation signal (HP-^He) is obtained. 

W To ensure that sufficient magnetization would be left after all radiofrequency (RF) 

pulses used in the imaging pulse sequence, a nominal constant flip angle of 12° was applied. 
S The repetition time (TR) was selected to be as short as possible as to reduce depolarization due 
O 20 to the presence of O2. The echo-time (TE) was also chosen to be as short as possible in order to 
□ minimize unwanted echo attenuation due to the high diffusion of the noble gas. The imaging 
''^ parameters were as follows: TR =18 minutes, TE = 3 minutes, matrix size 256 x 128, field-of- 
view (FOV) = 40 x 40 cm, shoe thickness = 7 mm, and number of sUces = 3. The voxel 
resolution was 1.5 x 3 x 7 mm. 
25 AST Perfusion Imaging . The AST perfusion data were acquired using a technique that 

had been previously described by Roberts et al, Rad 212(3):890-895 (1999). Specifically, a 
modified three-dimensional, compact gradient-echo pulse sequence was used used to acquire 
perfusion image data showing the radiofrequency (RF) and gradient (Gx, Gy, and Gz) 
waveforms (FIG. 1). The following parameters were used: TE = 716 |lis, TR = 50 ms, 7 mm to 
30 10 mm sUce thickness, 40 cm to 48 cm FOV, and 256 x 64 x 16 matrix. The voxel resolution 
was 2x8x8 mm. A firactional echo was acquired using a 64 kHz receiver bandwidth. 
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An arterial labeling RF pulse was applied in the presence of a gradient to provide 
continuous magnetic labeling of pulmonary arterial blood flow. Trapezoidal gradient and RF 
pulses were placed in the interpulse interval in order to label blood flow according to the 
principle of flow-induced adiabatic fast passage (Dixon et al, Magn. Reson. Med. 3:454-462 
5 (1986)). The offset frequency of the labeling pulse was adjusted in each case so that the 
location of the plane of arterial labeling intersects a pulmonary artery. The pulse, having an 
amplitude (Bi) of 24 mgauss and a duration (Ti) of 40 mseconds corresponding to a duty cycle 
of 80%, was applied during the time between the crusher/refocusing gradient pulses (C) and the 
beginning ofthe next excitation. The flip angle of the excitation pulse, a, was 30"^. The time 

10 axis is not to scale. 

The pulse sequence was triggered to the respiratory cycle to eliminate artifacts from 
breathing motion. The arterial labeling RF pulses and gradients were applied continuously in 
order to maintain a steady-state. Data acquisition was performed only during the end- 
expiratory phase ofthe cycle. Sixteen slab-encodings were acquired at end-expiration, 

15 requiring a temporal window of 0.8 seconds. Magnitude image data were acquired in the 

presence of arterial labeling ("inversion"), and under control conditions in which the amplitude 
ofthe label RF was set to zero ("control"). The acquisition ofthe entire perfusion data required 
less than 10 minutes per lung while the subject was breathing normally, but it could have been 
repeated without adverse effect on the data. 

20 In the normal subjects, the ventilation and perfusion signal was homogeneous 

throughout the lung fields as shown in FIGs. 2A-2C. However, the gas distribution and 
perfusion signals in patients with emphysematous lungs were markedly non-imiform (FIG. 3), 
and paralleled radiographic and nuclear V/Q scan abnormalities. The hyperpolarized He gas 
distribution and perfusion signals appeared relatively uniform throughout the lung allografts, 

25 with apparent defects caused by the heart and large vessels. Sagittal perfusion images using the 
arterial spin tagging (AST) method showed bright signals, corresponding to increased perfusion 
in the lung allografts. The emphysematous right lungs revealed poor perftision signals that 
were patchy and heterogeneous. A high level of ventilation and perfusion matching, confirmed 
with histogram analysis, was maintained in the single lung transplant recipients. 

30 A lower lobe wedge shaped perfusion defect, shown in FIG. 4, was observed in the 

patient with a large puhnonary embolus (PE), which was documented by computer tomography 
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imaging (CT angiography) imaging in the right lower lobe pulmonary artery. Although 
treatment with anticoagulants had been initiated when the PE was diagnosed several hours prior 
to imaging by the present metiiods, there remained a marked decreased perfiision signal in the 
lower lobe (shown witli arrows in FIG. 4). 

5 In sum, the present invention demonstrates that it is feasible to combine the use of a 

hyperpolarized noble gas (for ventilation) and arterial spin-tagging (for perfusion) to produce 
high-resolution, three-dimensional images of pulmonary V/Q in human subjects using a 
noninvasive and nonradioactive MR technique. 

Each and every patent, patent application and publication that is cited in the foregoing 

10 specification is herein incorporated by reference in its entirety. 

While the foregoing specification has been described with regard to certain preferred 
embodiments, and many details have been set forth for the purpose of illustration, it will be 
apparent to those skilled in the art without departing firom the spirit and scope of the invention, 
that the invention may be subject to various modifications and additional embodiments, and 

15 that certain of the details described herein can be varied considerably without departing from 
the basic principles of the invention. Such modifications and additional embodiments are also 
intended to fall within the scope of the appended claims. 
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